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Abstract 



We present a Thomas-Fermi treatment of resonant incoherent scattering 
of low-intensity hght by a dilute spatially confined Bose-Einstein condensate. 
The description gives simple analytical results and allows scattering data from 
finite-size condensates to be interpreted in terms of the properties of the ho- 
mogeneous BEC-system. As an example, we show how the energy dispersion 
of the elementary excitations can be measured from scattering by a finite-size 
atomic-trap condensate. As a second example, we point out that a near- 
resonant scattering experiment can observe quasi-particle creation caused by 
particle annihilation. 
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Typeset using REVT^ 



For most applications of the atomic-trap Bose- Einstein condensates [|T|] - it will be 
beneficial to have a highly populated condensate, and recent experiments have succeeded in 
increasing the number of condensed atoms 0]. It is fortunate then, that in this limit (for 
bosons interacting through an inter-particle potential of positive scattering length) many 
quantities can be calculated analytically. The simplifying assumption is the Thomas-Fermi 
approximation Jl^, which presumes that the local behavior of the BEC is similar to 



that of a homogeneous BEC with a chemical potential equal to the local effective chemical 
potential /i(r) = jiT — ^(i"), where [it is the chemical potential of the trap system and V{r:) 
the trapping potential. 

In this letter, we describe incoherent light scattering in a Thomas-Fermi approximation. 
The result, with a proper understanding of its limits, offers insight that will be useful in 
interpreting experimental scattering spectra. The formalism describes the scattering event 
as scattering of light by quasi-particles and includes the effects of recoil and intermediate 



state energy neglected in the off-resonant limit reported by Javaneinen ||T6|. 

The dynamical Thomas-Fermi approximation of this paper consists of describing the 
fluctuations that cause an incoherent scattering event to occur near a position r, by the cor- 
responding fluctuations in a homogeneous system of chemical potential /i(r). In this picture, 
the differential cross section for incoherent scattering, d'^ainc/dQdu, where dQ represents the 
infinitesimal solid angle and du the infinitesimal energy range {h = 1, in our units) over 
which the scattered particles are detected, reduces to an integral over a cross section density 



d'^a- 

^ "inc ^ I 

dVtdoj 



(P(^mc,hom/dflduJ 



(1) 

^l=^J.{r) 



V 

where d'^ai^chom/d^lduj is the cross section for incoherent scattering from a macroscopic, 
homogeneous system of volume V and chemical potential /i(r). 

We use a second quantized representation in which Ck, cj^ denote the annihilation and 
creation of a single-atom state with the atom in the atomic ground state (the trapping state) 
and its center-of-mass in the plane wave state of wave vector k; Ck, cj^. denote the correspond- 



ing annihilation and creation operators for atoms in the excited (resonant) atomic state. In 
the resonant scattering process, the interaction with the electric field of the incident light, 
Ein hn exp(z [kj„ ■ r — cUmt]), where Ein represents the intensity and the polarization of 
the incident light, promotes an atom to its excited atomic state. The excited atom subse- 
quently deexcites, creating a photon of electric field Eq iout exp(— i [kout ■ i" — ^outt]), where 



Eq is the single-photon intensity (in Gaussian units, Eq = y 27rA;c/V, where k is the resonant 
photon wave number). The corresponding absorption and emission Hamiltonian operators 
in the interaction picture are : 

^(abs)(^) = -Ein (d ■ iin) exp(-z^i„t) ^4+k,„(^)ck(t) , 

k 

H(-^){t) = -Eo (d ■ eoutT expizuoutt) E cl.^.Jt)d^it), (2) 

k' 

where d denotes the resonant atomic dipole moment. 

The amplitude of a state \F) of the scattering system, a time T after the system was in 
its initial state \in), is, to lowest order in perturbation theory, equal to : 



/■T/2 i-T/2+t 

A{\in), \F)-T) = - / dt dT{F\H^'''"\t)H^^'''\t - r)\tn) , (3) 

J-T/2 Jo 



where the photon absorption occurs a time r earlier than the photon emission ||18[. In the 
limit of low-intensity incident light, each photon is scattered by a system of ground state 
atoms. Indicating the many-body vacuum of excited-state atoms, |0), in the initial and final 

states, \in) \in) |0), \F) \F) |0), we obtain 



A{\in), \F); T) = -E,^ Eo (d ■ e,„)(d ■ iout)* dt dr 

J-T/2 Jo 

E (^l4'-k(^)ck(t - r)|m) (0|cu'(t)4+k,„(i " ^)|0) exp(^ [ujoutt - a;,„(t - r)]) . (4) 



k,k' 



The vacuum expectation value in (^ is equal to (0|ck'(t)4+,^. (t — t)|0) = 5k',k+ki„ 



X exp —i [iUo + -Ek+ki„ — i'y/2j T , where we write the energy of the excited atom as the 
sum of the atomic excitation energy, c^o, and the kinetic energy, ii^k+ki„, and where 7 is 
the width of the excited atomic state, 7 = {4/3)k^d'^. The upper limit of the r-integration 
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interval, may be replaced by oo if T/2+t ^7 ^. Since we shall take the limit T 
change the r-interval accordingly, and we introduce a scattering operator Sq{t): 



rT/2 

A{\in), \F)- T) = -Ein Eo (d ■ Q„)(d ■ iout)* / dt exp(zcut) {F\S^{t)\in) , 

J-T/2 

where S^{t) =Y. 4+q(^)ck(^ - r) exp(-z \ujq + ^k+k,„ - - «7/2 



00, we 
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where q is the momentum transfer, q = kj„, — hout, and u the energy transfer u = ujin — ujout- 
The scattering rate is the ratio of the square of the amplitude @ over T, summed over 
all |F)-states and thermally averaged over the |m)-states, in the limit of T ^ 00. The 
differential cross-section d'^ahom/dQdu is the product of the resulting rate with the ratio of 
the scattered particle final state density, pout, over the incident particle flux, Jj„. For photon 
scattering, pout = [^/(2vr)^] x {k'^/c), and Jin = (Ef^/E^) x (c/V), so that 



d'^aiio^/dQdu = lim — 
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where d = d/d, and the ( )-brackets denote the thermal average over the (m| |m)-matrix 
elements. 

To describe scattering from a dilute BEC of low depletion, (A^ — No)/N << 1, where 
is the total number of atoms and Nq the number of condensed atoms, we work in the 
Bogoliubov approximation. In this scheme, we treat the zero-momentum operators as c- 
numbers, replacing them by -\/iVo, and we keep only terms proportional to \/Nq: 



S^it) 
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dr cUt) exp(-z cjq + ^k,„ - i^in - ^7/2 
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The first term on the right-hand side of (|^ represents the scattering event in which an atom 
is taken out of the condensate, and the second term represents the event in which an atom 
is put into the condensate, each process transferring a momentum q to the many-boson 
system. Since kj„ — q = kout, the excited atom momentum equals the resonant photon wave 



number k in both processes, and we absorb the recoil energy Ek into the definition of the 
detuning: A = — ojq — Ek- 

Finally, we perform the Bogoliubov transformation to quasi-particle operators 6, 6^ : 

c\{t) = cosh(crq) - sinh((Tq) 6-q(t) , 

c-dit -t) = cosh((Tg) 6_q(t - r) - sinh((Tq) 6^(t - r) , (8) 

where the canonical nature of the boson b-operators is ensured by writing the coherence 
factors as cosher and sinha. With the time- dependence of the creation and annihilation 
operators, &_q(t — r) = &-q(t) exp{iEgT), bl^{t — r) = fojj(t) exp(— ii?gr), we obtain for (0) : 

cosh (7g sinh ag 

^ [A + Eg + i-f/2 ~ A + i-f/2 

which describes the response of the BEC to the scattering as creating and annihilating quasi- 
particles. Each term of (P) represents a different process that creates or annihilates a quasi- 
particle. For example, the scattering can create a quasi-particle (6)^) either by scattering an 
atom from the condensate into a state of momentum q (c^) , giving an amplitude cosh cfg for 
quasiparticle creation or by scattering an atom from a state of momentum — q (c_q) 
into the condensate, giving an amplitude sinhcr^ for quasi-particle creation. The resonant 
denominators are the usual (complex) energy-differences between the initial and intermediate 
state: in the first process, the energy-difference is ujin — {ujq + E^ — i'^ (2) = A + i'j/2, whereas 
in the second process, the intermediate state contains an extra quasi-particle and the energy 
difference is equal to — (cuq + Ek + Eg — /2) = A — Eg + i'y /2. 

At zero temperature, (6Lq(t)6-q(0)) = 0, and (&q(i)&q(0)) = exp{—iEgt), so that 

^ f dtexp{iujt){Sl{t)Sd{0)) ^ No 6{uj - Eg) 

giving a spectrum that consists of a single peak. The scattering processes neglected in the 
Bogoliubov approximation, give rise to an additional background in the spectrum. The ratio 
of the peak intensity to the integrated background intensity is roughly proportional to the 
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A + i^/2 A-Eg + t-f/2 



(10) 



depletion. The finite-temperature generalization of (|TOD can be achieved in a straightforward 
manner by introducing temperature-dependent quasi-particle occupation numbers z/g, so that 
(6^q(t)6_q(0)) = Uqexp^iEqt) ctc... (scc for example ref. ||T9[] ). 

The value of is determined by minimizing the free energy. At T=0, the Bogoliubov 
approximation gives tanh2(Tg = /i/ (g^/2m+/i), where /i is the chemical potential which, for a 
condensate of atoms interacting through a potential of scattering length a, is yU = {An /m)apo 
where po is the condensate density. Furthermore, Eg = {q^ /2m + pY — f^'^, and with 
and (|10|) , we obtain the following expression for the cross section density in (|l]): 



'^^O"inc,hom/ dVtduj 

V 
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{q^/2m + p + Eg)/Eg 



{qy2m + p-Eg)/Eg 



A + i^/2 A-Eg + i^/2 ' 

where Eg and p are r-dependent. The energy-conservation factor, 6 {u ~ Eg{r)) implies that 
an energy interval {uj,uj + duj) probes the condensate region 00 < Eg{r) < 00 + duj. In a 
magnetic trap, the excited atom experiences a potential energy due to its magnetic moment. 
We account for this effect by making the detuning position-dependent: A(r) = A — aV{r), 
where a is the ratio of the excited atom potential and the trapping potential. Finally, for 
the sake of simplicity, we specialize to a spherically symmetric harmonic oscillator trap, 
V{R) = {ujt /2){R/ L)'^ , where is the trapping frequency and L the extent of its single- 
particle ground-state, L = l/y/mur- In the Thomas- Fermi approximation, the condensate 
density is ( fTl, [|ll) po{R) = [pt/^] x [1 - (R/Ro)'^] 0{R-Rq), where Rq is the condensate 
radius, Rq = L{lhaN / L)^^^ and pr = {ujt/2)(Ro/L)'^. In performing the spatial integration 
(H), spherical symmetry reduces the expression to an integral over the radial distance R, and 
we substitute R by the effective chemical potential p, R = Rq ^Jl — p/ px- With this substi- 
tution, po{R) ^ /i/^, A(i?) ^A-a[pT-p] and 5{u - Eg) ^ 6 {p - Pg{uj)) \dEg/dp\-\ 
where Pg{uj) is the effective chemical potential at the positions where Eg is equal to uj, 

/^g(^) = I ^272^ - 2^ ' and ^ = 

Pg{uj) 



^-^r^, we find 
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The single peak at cu = Eg in the spectrum of the homogeneous BEC ([TT| ) is broadened 
to a feature from uj = g^/2m to = 



{q^ /2m + /iy)^ — /^T- The region in parameter- 
space where the difference in intermediate state energies can be neglected, either because 
|ci;| ^ |A — a{iiT — oi' 1^1 ^ 7/2 is the 'fast collision' regime. As in the off- 

resonant limit [0 |T^, which is part of the fast-collision regime, fast-collision resonant 
scattering data contain the same information as non-resonant scattering data |T^, giving a 
cross-section proportional the dynamical structure factor of the scattering system: 



(in the fast collision regime 



where StfIq, t^) 



dV 



(37/4k)(d-ein)(d-eout) 



STF(q,tu) 



2{jJt^ V aL^ 
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(13) 



is the dynamical structure factor of the condensate, calculated in the above Thomas-Fermi 
approximation. In Figure 1, we compare the cross-section (|T2|) and the fast-collision approx- 
imation (^) for three different values of the detuning. 

We also point out that the dependence on the energies of the intermediate states implies 
the interesting possibility of observing the creation of a quasi-particle caused by particle 
annihilation. The A — + i7/2-denominator corresponds to an intermediate state of higher 
energy than the A + i7/2-intermadiate state, indicating quasi-particle creation^ whereas the 
intermediate state was formed by the removal or annihilation of a boson (by exciting it to 
a different atomic state). The observation of the difference in energy-denominators can be 
accomplished by varying A while keeping uj constant, and requires near- resonant detuning 
which, unfortunately, can make the condensate optically thick. Nevertheless, experimental 
techniques such as resonating on different atomic transitions and/or using the polarization 
of the photons (scattered by spatially oriented dipole moments), can reduce the optical 



thickness. 

Finally, to understand in what sense the scattering data can be interpreted in a Thomas- 
Fermi manner, we need to investigate its limitations. A necessary condition for the validity 
of the dynamical Thomas-Fermi description, is the validity of the static Thomas-Fermi 
description. At T = 0, in a harmonic trap, a Thomas- Fermi condensate satisfies Rq » L, 
or equivalently, ^ » Hujt ( [|n|], ||l2l, |^). Furthermore, the inability of our dynamical 



Thomas-Fermi approach to describe the discrete spectrum of a finite-size BEC, indicates 
that even for a Thomas-Fermi condensate, the dynamical Thomas-Fermi results should be 
interpreted within certain limits. For the scattering problem, we need to realize that only 
fluctuations confined to limited regions in space and time are well-described by a Thomas- 
Fermi approach. The first constraint pertains to the excited-atom propagation, which is 
only well described by means of a position dependent detuning if, during its lifetime, the 
excited atom experiences a change in potential small compared to /17, ([|F ■ v| /•y] << ^7) 
where F is the force and v the velocity of the excited atom. A second condition pertains 
to the Thomas-Fermi description of the BEC-fiuctuations. Consequently, we require that 
q > where ly indicates the scale on which /i(-R) varies (for the harmonic oscillator trap 
we can choose ly = Rq/S ; fiesiR) varies by approximately 10 % from _R = to -R = Rq/S). 
This condition is somewhat stricter than the condition to observe incoherent - as opposed 
to coherent - scattering, q > Rq^. Nevertheless, with a chemical potential /i^ ~ 100(u;t'/2) 
(realistic in present-day BEC technology) Rq ~ lOL and — 10/L, where is the inverse 
of the local coherence length in the middle of the trap k^ = 2^JJlTm. We find then that the 
condition q > ^ 3/i?o — 0.03/Cc; leaves almost all of the interesting part of the dispersion 
(expected to be phonon-like up to kc) to be explored. As for the temporal constraint, the 
'locally homogeneous- like' excitation picture breaks down on the time scale that it takes 
the BEC-response to be affected by its inhomogeneity. Assuming that an excitation in the 



middle of the trap propagates at the local speed of sound, c = \J ^t/^ = 'x -Ro/v^, we 
can estimate the relevant time scale as ty ~ ^'"^^"^ . Reducing the frequency resolution in 
the scattering spectrum to Acj ~ ~ ciJt x (Ro/h), restricts the scattering probe to short- 
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time {t < ty) 'homogeneous-like' fluctuations. Thus, the Thomas-Fermi scattering spectrum 
should be interpreted as a 'smooth' version of the real spectrum and we should compare 
intensities integrated over frequency intervals larger than or equal to with experimental 
results. These estimates were made in the middle of the trap where the Thomas-Fermi 
description works best, and which is probed on the high-frequency side of the cross-section, 
whereas the low frequency-region, uj ~ /2m, probes the edge of the condensate where the 
Thomas- Fermi results cannot be trusted. 
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Figure 1. Differential cross section for resonant light scattering from a trap with the 
parameters shown in the figure caption. The plots show the cross section as a function of 
energy transfer, at fixed momentum transfer q = ^muT, for three values of the detuning 
A. The full line shows the Thomas- Fermi calculation and the dotted fine shows the fast 
collision approximation result. As the detuning increases, the fast collision approximation 
becomes more and more accurate. 
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